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I. Objectives of the Manual- 

The orbiting Geophysical observatory (OGO) satellites move 

in the ionosphere and the magnetosphere of the earth and carry on 

joI ; -rot- 

board various experiments to obtain scientific informations on the 
geophysical properties of the atmosphere pertaining to the solar- 
terrestrial relationships. Although the satellites are designed to have 
booms to carry detectors which extend away from the gasdynamical 
disturbances generated by the moving satellite, it is not always 
possible to achieve the goal in view of the various orientations, with 
respect to the direction of motion, that a satellite must assume during 
its orbital motion. This is further complicated by the fact that the 
extent of the disturbed region for a tenuous plasma such as the ionospher 
is quite large. It is, therefore, almost unavoidable that a detector 
. or intake orifice happens to situate in a disturbed flow field which is 
generated by the other parts of the moving satellite. Consequently, 
it becomes imperative to be aware and understand the disturbed-flow 
effect on the measurements in order to minimize the error of 
measurement so incurred. 

In determining the field of disturbance around an orbiting OGO, 
one must make considerable flow analysis for both the neutral and 
the charged particles of the ambient atmosphere pertaining to, various 
complex geometrical configurations of an OGO. The task can be 
exasperating and even the presentation of the calculated results of the 
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flow field becomes intractable as the guidance for its users. To 
circum\ ent such difficulty, a practical manual for the users is devised 
to provide the basic informations of flow fields which can be synthesized 
for a given geometrical configuration with minimum effort. A manual 
must satisfy two basic requirements: (1) the informations on field 
variables, e.g., particle densities, electric field potentials, etc., 
must be readily available from tables and charts without undue 
computations on the part of the users; (2) the collections of the basic 
informations on the field variables of interest must be versatile 
such that it is possible to synthesize for various geometrical config- 
urations and orientations of the orbiting OGO satellites. 

In view of the basic characteristics of restricted zone of action 
for supersonic flows, such as those of the orbiting satellites, it is 
found most fruitful to achieve the above-mentioned objectives by 
presenting the field variables in charts under various ambient conditions 
for simple geometrical bodies, e.g. sphere, cylinder, disks, etc., 
of which OGO satellites are composed. The superposition of appro- 
priate basic flow fields with due account of the mutual interference 
from the proximity of the neighboring elements can be made to estimate 
the flow field of interest, 

II, Preliminary Informations 

• The qualitative and quantitative nature of the disturbance 
induced by an orbiting satellite must depend, among other things, 
upon the ambient atmosphere in which it is moving and the body 
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parameters characterizing the surface potential, size, and speed, 
etc. , of-the moving object in question. Two tables are provided, 
one of which describes the averaged properties of the upper ionosphere; 
the other, the properties of typical satellites. 
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Distance from . 
surface of the earth 

Concentration of 
neutral particles 
(no. / cm 3 ) 

Concentration of 
electrons and ions 
(no. / cm 3 ) 

Ambi ent t emp . 

(°K) 

Mean molecular wts. 
of the particles 

Effective Coll, 
frequency 
(electrons, ions) 

Mean free path 
of neutral particles 
(cm) 

Mean free path 
of an electron among 
ions (cm) 


Thermal speed of 
neutral particles and 
ions (cm/ sec) 

Thermal speed of 
electrons 

(cm/ sec) 

Larmor radius of 
electrons (cm) 

Larmor radius of 
ions (cm) 


Debye Length (cm) 


2'50'-400 Km 

600-700 Km 

— 3000 Km 

.1.0 10 - 10 9 • 

10 7 


5 x 10 5 - 5 x 10 6 

10 5 

10 3 

— 1000 

— 2000 

-4000 

— 24 

— 1 6 

— 2 

(1 - 3) x 10 3 

(1 - 3) x 10 2 

10 

(1 - 7) x 10 s 

(1 - 5) x 10 7 

2 x 10 14 

— 10 4 

~10 5 

— 10 6 

— 10 s 

— 1.5 x 10 5 

— ' 6 x 10 5 

~2 x 10 7 

— 3 x 10 7 

-4 x 10 7 

-v 3 

— 4 

-10 

— 5 x 10 2 

-8 x 10 2 

-10 3 

— 0.5 

"'■'0.5 

— 3 
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Table II a - Propertie 

s of Satellites 



Distance from 
surface of the 
earth 

250-40.0 Km. . 

600-700 Km 

3000 Km 

Body velocity 
(cm/ sec) 

O 

i—H 

i, ' = ' 
1 

CO 

(8 - 11) x 10 5 

(8- 

10) x 10 5 

Speed ratio, 
body/ thermal(ion) 

8 ■ - T’l 

5 - 7 


1 

Table II b 

- Accommodation of Particles on a Solid Wall* 


Neutral particles 
accommodated 

Air He 

A He 

A 

He 

Surface 

polished Tungsten 
Al. wire 

flashed 

T , wire T . wire 
flashed not 

flashed 

T. wire 
not 

flashed 

glass 

or 

0.87-0.95 0,17. 

0.82 0.53 

1.00 

0.31 


*S. Flugge(ed) Handbuch der Phys. 8/pt,2 (Mech. of Rarefied Gases), 1963. 
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III. Neutral Particle Distributions near a Moving Satellite 

It can be seen from Table I on the properties of the ionosphere 

that the mean free path of the ambient particles in question are many 

orders larger than the characteristic dimension of the satellite. 

Consequently, the neutral particles, e. g. molecules and neutral atoms, 

of the atmosphere interact with the satellite according to the theory of 

the free molecular flow. Reflection of the incident stream to the 

satellite creates an excess of particles in front of the body where an 

enhancement region is formed. Behind the body, a rarefaction region 

develops since the rapidly moving body sweeps the particles, the 

ambient particles cannot fill the cavity aft the body immediately. The 

length of the rarefaction wake increases with the velocity (V ) of the 

o 

body. An important parameter here is the speed ratio, V /v where 

o th 

V th is the thermal speed of the ambient molecules that fill the wake. 

Note that the thermal random motions of the molecules tend to damp 
out the disturbances generated by the motion of the body. 

To determine the flow field around a. body which can be decomposed 
geometrically into simple basic elements of sphere, cylinder and 
disk, one may use a principle of flow synthesis wherein the flow 
field of the combined configuration is obtained by patching the component 
fields and taking the mean values in the overlapped regions. 

Ilia. Neutral Particle Distribution in Front of a Moving Body 

The concentration (or number density) of the particles in front 
of a rapidly moving body such as a satellite are contributed from 
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two sources: (1) the undisturbed ambient particles of concentration 
n 0 prior to their impacts with the body, (2) the particles reflected 
from the body upon collisions. The density and shape of the cloud 
of reflected molecules may be of importance in many OGO- experiments . 
The number density of the reflected particles in this cloud decreases 
as the square of the distance from the body, due to the essentially 
spherical nature of the distribution of the- reflected particles with 
respect to the body. The concept of free molecular flow will be 
applied to the analysis of this cloud. 

It is noted that the quantitative description of the cloud of 

reflected molecules depends critically on the nature of the reflection 

mechanism at the surface of the body. Although the reflection mechanism 

for very high speed molecules may not be adequately described by 

the classical models of specular and diffuse reflection, it is felt, 

however, for the mean molecular energy of impact on OGO- satellites 

close approximations are expected with the use of these models and 

judicious choice of accommodation, coefficient (a) which is defined 

as the fraction of the mean translational energy of the incident 

molecules which is transferred to the surface. The quantity or is 

thus a measure of the degree to which the molecules have their mean 

energy accommodated to what it would be if the molecules were 

reflected with a Maxwellian distribution corresponding to the surface 

temperature (T ). For perfect accommodation a = 1 also called 
w 

diffuse reflection; while for the case of vanishing energy exchange 
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with the surface) a; = 0 , it is called specular reflection where the 
molecules reflect as a light beam on a mirror. A statistical 
approximation of the true reflection which obviously depends on the 
molecular energy of impact upon the surface and the smoothness of 
the surface finish among other things. Typical values of a are shown 
in Table lib. The approximation for the density of reflected 
molecules can be made with n = cyn(diffuse) + (1 - cr)n(sp). 
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(i) Reflection from a Small Surface Element 


Consider a stationary surface element making an angle 0 with 

the free stream velocity (V ). The free stream density (or ambient 

density) and temperature are n^ and / respectively. The 

temperature of the surface element is T . Let v, = (2kT/m.) 1/<2 

w th ' 

and S = Vq/v^ where k = Boltzmann constant and m = molecular 
mass. For diffuse reflection {a = 1), the number density of reflected 
molecules at the wall is 


j _i 

n w ~ ~2 n 0 ^ T 0 ^ T w^ 2 {? xp (“S 2 sin 0) + \[tt S sin 0 (1 + erf [S sin 0]^) 


x 


, , 2 r~ -u 2 du 

where erf x = — p I e 
Vir J 
o 

It is noted that in OGO flight, S is in the order of 10, unless sin 0 
is very small, the following approximation is valid 

n w~ V o Sme. 

When the molecular specular reflection at the surface is assumed, 
a - 0 , the angle of reflection is equal to the angle of incidence for 
any incident molecule. Since available evidence suggests that specular 
reflection is likely to be of importance only in very high' speed flows 
(S » 1), the analysis will be restricted to this case. The number 
density of the reflected molecules at the surface is very nearly equal 
to (for 0 > 0) 


n = n . 
w o 
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(ii) Reflection from a Flat Plate 

Consider a stationary flat plate with, the free stream velocity 

(V^) normal to the plate (Figures 1 and 2). The shape of the cloud 

of the reflected molecules as indicated by the contour of constant 

density n = 10 4 n for either the case of nearly diffuse reflection 
wo 

(Figure 1) or specular reflection (Figure 2) under different speed 
ratios (S) are shown. 

(iii) Reflection from a Cone at Zero Incidence 

Consider a stationary cone with a base radius R and a semi- 

o 

vertex angle (3 exposed to a stream of velocity (V ) at zero incidence 
(see Figure 3). With diffuse reflection the number density of the 
reflected molecules at a point P (r, cj>), where r, <j> are the polar 
coordinates with origin at the tip of the cone, is given by 

= w r\ v./h 7 , 

n r 4 it sin R x J y T 1 

w 

where 

x = exp (-S 2 sin (3) + \Trr S sin (3 [l + erf (S sin (3)] 

and 

Fj = cos tb sin (3 - cos 1 (cot tan j3)J- + sin <j> cos (3 sin ^rr - cos' 1 (cot <p tan 

Noted that S is large (about 10) for OGO, hence, for a cone 
with sin (3 > 0.2 an approximation for n can be used 



With specular reflection; the number density of the reflected 
molecules at P(r, <j>) becomes 
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n S 
o 


n 

r 



(sin<j> sin 2(3) 2 exp£-2S 2 [l 


+ cos ( 2(3 + <j>)]} 


x 


erf [S (sin (j> sin 2(3) 2 ^ - cos " 1 (cot<j> tan ( 3 )^] 

Figure 3 presents some results of reflected molecules for both 
diffuse and specular reflections at the cone. 

(iv) Reflection from a Sphere 

Consider a sphere of radius R^ which is exposed to a stream 
of velocity (V^) as shown in Figure 4. The expression for number 
density (n^) of the reflected molecules at a point P(r, <J>) can no 
longer be given in closed form. With the assumption of S » 1 
which is valid for OGO, the number density of reflected molecules 
with an assumption of diffuse reflection can be reduced to a simple 
integral which must be integrated numerically: 


n /R \ 

■fCf) J 


n 


2 2L 
•2 


sin 2 £ F 2 (<f>, 1) d£ 


where 


- cos 4 cos £ ^jr - co s ' 1 (cot <j> cot + sin 9 sin £ sin £tt - cos x (cot 9 cot £ 
With specular reflection, the number density of the reflected 
molecules at P(r, 4 ) is 

n r = WitCr) J 2 exp {“ 2S2[1 ' 003(26 - 


erf [S (sin 9 sin 2 £) 2 - cos ' 1 (cot 4 . cot £)} ] d£ 

which also must be integrated numerically. 
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Some results of the numerical integration for the number densities 
of reflected molecules for both diffuse and specular reflection at 
a sphere are shown in Figure 4. 

(v) Reflection from a Circular Cylinder 

Consider a cylinder of radius R q and length h having its 
axis aligned with a stream of velocity (V ) . The number density of 
the diffusely reflected molecules from the curved surface at a point 


P(r 5 <j>) 


n 

r 



sin <j) 


Illb. Neutral Particle Distributions Behind a Moving Body 

The concentration (or number density) of neutral particles 
behind a moving body is also attributed to two sources; (1) those 
filling into the cavity behind the rapidly moving body; (2) those 
reflected from the surface of the body. It is noted, however, that 
the second source of input decreases as the speed of the body increases. 
For amoving body, such as OGO, the contribution of the reflected 
molecules in the wake becomes insignificant except in special situations 
such as that discussed in section IV of this manual. 

(i) Neutral Particle Density in the Wake Behind a Sphere or a Circular Disk 
As a sphere moves rapidly in a tenuous atmosphere, it 
continuously bores out a cylindrical cavity in the medium. The 
neighboring ambient molecules must move in to fill the cavity which 
remains, however, underpopulated. Consider a stationary sphere 
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of radius R^, the center of which is at the origin of the cylindrical 

coordinate system, (z, p , 9) with z aligned with the free stream 

velocity (V^), the filling action of the ambient molecules is provided 

by their random thermal velocities.. The wake structure is obviously 

axisymmetrical, hence, the position-dependency of number density 

in the wake is two-fold (z, p). To simplify the expression, the 

position coordinates are measured in units of R , the radius of 

o 

the sphere. 

The expression for the number density in the wake of a moving 
sphere must be specified differently for different regions of the 
wake. Where p < 1 , the number density is given by the following 
integral which must be integrated numerically: 


~ (1 - erf S) + (4-n-r 1 / 4(1 + A 2 f 2 exp 

° o l 


1 + A 2 


x 


1 + erf 


AS 


yr+ a 2 _ 


d d 


o 


where 


A - z(l - p 2 ) 1 [( 1 - p 2 sin 2 <J>) 2 - p cos <j>] 


2 r -u' 


erf S - —h J e du 


while for p > 1 
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n(z, p ) 

n 

o 


arc sin (1/ p ) x 

o2 

1 + (4-ii-r 1 j {A(l + A 2 ) -2 exp 

-S 

1 + A 2 

-arc sin (l/p) 


1 4-erf 



i 

B(1 + B 2 )~ 2 



B‘ 


1 +B 2 


■] 


1 + erf 


bS 


ri 


+ B" 



where 


B = - z(l - p 2 ) 1 [(1 - p 2 sin 2 <p) 2 + p cos <j>] 

It is noted that in the above analysis, the contributions to the 
cone entr ation of particles in the wake due to those reflected from, 
the surface of the moving body have been ignored in view of 'S » ] 
for OGO. Consequently, the result is independent of the shape of ’ 
the body except the contour of its maximum cross section facing 
the stream. Therefore, the result is identical for either a sphere 
or circular disk in motion except in the semi -spherical region 
r — B. „ Some results of the numerical computation of number 
densities of neutral particles in the wake behind a sphere (or circular 
disk) are shown in Figures (5), (6), and (7). 

(ii) Neutral Particle Density in the Wake of a.,Cylinder or a Rectangular Disi 

Consider a cylinder of radius R and length 2R with its 

o x 

axis on the x-coordinate and the free stream velocity (V ) along 
the z- coordinate. Let the position coordinates be measured in 
units of R^ ; velocities in units of v = V2RT/m . The number 
density at a point (x, y, z) in the wake of the cylinder is 

- - 00 T y 

“ = 2 J exp[ - (v --S ) 2 / I - - [ erfM - erf N][erf P - erf Q]l dv 

o -oo Z 1 4 j z 
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where 


M = (v /z )(x - R ) N = (v /z)(x + R ) 

z x z x 

P =• (v / z)(y - R ) Q = (v / z)(y + R ) . 

z o z o 

As in the case (i), the contribution due to molecules reflected 

from the surface of the cylinder has been ignored in view of the 

high speed ratio (S » 1) of OGO. Consequently, the above 

expression is also valid for a rectangular disk with the free stream . 

(V ) normal to the surface of the disk 2R x2R .. 
o ox 

Some results of the neutral particle distributions in the wake 
behind a cylinder (or disk) are shown in Figures (8), (9), (10), (11), 
(12), (13), (14), (15), and (16). 
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figure is - Concentrations along 
radial directions at 0°, 27° 

45 °, 60° with V 0 behind a 
cylinder at S =8. 








V. 


Charged Particle and Electric Potential Distributions Near a Moving 
Satellite 

The problem, of charged particle distribution in the vicinity of 
a moving body is much more complicated than that of the neutral - 
particles. The fact that the particles are electrically -charged and 
the satellite in motion is a conducting body means that the latter 
will become electrically charged even if it is started with zero 
surface potential. In the near earth environment where the charge 
input to the satellite by the ambient plasma dominates over that by 
the photoelectric effect due to solar radiation, the -equilibrium 
surface potential is negative in the order of a few volts. This 
condition has been assumed in the following discussions of the 
charged particle distribution near a satellite. Such a condition 
could become invalid under special circumstances of the OGO 
flight. This possibility will be studied in a future program. 

One of the charged-particle effects is that the particle motion 
is coupled to the electric potential field which, iri turn, depends on 
the space charge contributed by the moving particles-. A solution 
for the particle distribution in the vicinity of a moving satellite 
must be self-consistent with the solution for the electric potential 
distribution of interest. The results presented in the following are 
self-consistent solutions. In the study of the plasma interaction 
pertaining to a satellite which moves with a mesothermal speed 
(much higher than the thermal speed of the ambient ions and yet 
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much lower than the thermal speed of the ambient electrons) it is 
convenient for study purpose to divide the region of plasma distur- 
bance into two parts: (1) the sheath which refers to the thin layer 
lying between the ambient plasma and the frontal surface of a 
rapidly moving body; (2) the wake which denotes the cloud of particles 
immediately behind a moving body. 

IVa, Plasma Sheath 


There are two important characteristics of the sheath associated 
with a negatively charged body moving at the mesothermal speeds 
that renders simplifications for its analysis: (1) cold-ion approxi- 

mation; (2) quasi-one dimensional nature of the sheath structure. 
These are recognized in the sheath study. 

A new concept using the Schrodinger wave equation together 
with the Poisson equation for the sheath of cold-ion plasma at meso- 
thermal speeds has been proposed and a solution for the sheath of 
an inclined charged surface of negative potential is obtained. Let 
the free stream velocity (V ) impinge on a large conducting plate 
of potential tp^ at an angle 8 to its normal which aligns with the 

x-axis; free stream ion density = electron density = n ; free stream 

o 

ion kinetic energy, E = — mV 2 , electron temperature, T . The 

z o e 

electric potential cjj(x) in the sheath is given 


$ 


x-X q = J (8Trn o )" 2 {2Ecos 2 0(Yl - e<|)/Ecos" z 7- 1) - kT [l - exp(ecj>/kT )] 


' 4 > 


where $ denotes the potential at x = x . Some results of the 
o o 


J Ul 
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electric potential and ion density in the sheath along an inclined 

plate are shown in Figures (17) and (18) wTaere S denotes the 

x 

component speed ratio. (For details see Appendix A)., 

The extension of the above analysis for a curved surface is 
feasible . It is noted, however, that the sheath behaves much like 
the boundary layer in fluid mechanics in that it varies much faster ' 
in the direction normal to the surface than along it (the quasi-one 
dimensional nature). This makes it possible. to segment the curved 
surface into short inclined planes for which the sheath theory is 
known. 

IVb. Distribution of Charged F.articles and Field Potential in the Near Wake 
With a satellite, e.g. OGO, moving at a mesothermal speed, 
the nature of the plasma disturbances of the wake differs considerably 
from that of the sheath in front of the body. When’ a body moves 
rapidly through a plasma, a momentary cavity is created behind 
the body. The ambient ions and electrons rush into the cavity by 
virtue of their thermal velocities. The ionic thermal velocity, 
though small compared with the free stream velocity, cannot be 
ignored in the wake analysis as done in the sheath analysis. 

■ It is the search for the self-consistent solution of the plasma 
disturbances of the near wake that much of our efforts have been 
directed at. 

To cope with the complex geometrical configuration of an 
OGO, we again consider the satellite as composed of simple geometrical 
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components and. take into account the interference effect due to the 

proximity of the neighboring members in the ensemble. The details 

of these new self-consistent solutions for the distributions of charged 

particles and electric potentials in the wakes of cylinders and spheres 

are described in the Appendices B and C. We shall present here 

the computed results in Figures (19 - 72) for spheres and cylinders 

under various parameters of satellite environments and ambient 

conditions. Each solution is prescribed by four parameters: 

(1) R/d, (2) T^/T., (3) V^/w, and (4) cb^. i. e„ The ratio of body 

size (R) and Debye length (d), the ratio of electron and ion temperature, 

the ratio of free stream velocity and thermal ion velocity, and the 

surface potential (normalized by kT./e), respectively. 

The common features of these solutions as shown in Figures 

(19 - 72) are as follows: (1) the mean free path of the ambient 

particles are many orders larger than the characteristic dimension 

(R) of the body which is, in turn, large compared with, the Debye 

shielding length (d =‘VkT / 4 tt e 2 n') of the. ambient plasma; (2) the 

velocity of the moving body (V ) is me so thermal, i.e. to say, it 

is much larger than the ionic thermal speed (v. = V2k T./m ' ) and 

i i i 

yet much smaller than the electron thermal velocity (v ="\/2kT /m ! )• 

e e e - 

(3) the surface potential is negative relative to the ambient plasma. 

The geomagnetic effect has been ignored. This is justifiable provided 

the ionic Larmor radius at the satellite altitude in question is large 

* 

compared with the body size (R), The electron density is assumed 
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to follow Maxwell- Boltzmann distribution, namely, n =n exp(T. a^/T 

e o 1 “ e 

where n denotes the ambient electron density: cb* , the local 
o 

electric potential (normalized); hence, n is not tabulated here . 

c 

It is of interest to describe some of the significant conclusions 
of the present calculations: (1) It is found that the distributions of both 
the particles (n. and n ) and the potential (<p=i<) are insensitive to 
the change of R/d as long as R/d » 1; e„g. the difference in the 
solutions for R/d = 20' and R/d = 100 is negligible. (2) As long as 
R/d 25> 1, there is always a pronounced potential well of negative 
voltage present in the wake at a distance of (1 - 2)R from the 
rear surface of the body. Even with R/d = 1, there is a potential 
well of mild intensity in the wake. The theoretical prediction of 
the presence of this potential well has led to the hypothesis of trapped 
particles in the wake which might explain the anomalous plasma 
oscillations in the wake. (For details see Appendix D.) (<j)*= e<j>/kT.) 

V. Illustration of Applications 

The procedure for the use of the present manual in prediciting 

the disturbances in distributions of particles and potential near a 

moving satellite can be illustrated by an example. Consider a satellite 

■ with a configuration as shown in Figure (72). It consists of two 

spheres A and C which are connected by a long cylindrical rod 

of radius R 3 , all of which has a normalized surface potential of 

d> = - 1 (kT./e) volts, 
o i 
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it is assumed mat me satellite 
is orbiting at an altitude of 600 Km. 

It is of interest to know the distributions 
of charged and neutral particles and 
also the electric potential near the 
satellite. 

The procedure to determine 

the quantities of interest can be 

stated as follows: 

(a) neutral particle distribution - 

The configuration of Figure (73) is 

resolved into three components: spheres A and C, and cylinder. B. 

It is noted from the figures of computed solutions of neutral particles 

Figures (5 - 16) that the parameter of primary interest is the free 

stream speed ratio V / v, . The thermal speed (v , ) of the ambient 

o th . r th 

particles is estimated, by using Table I, to be 1.5 x 10 5 cm/sec. - 

Thus, V /v,, can be easily obtained. Knowing V /v , , the 
o th ° o th 

particle distributions around the components A, B, and C can be 
determined by the use of formulas in Section II and Figures (5 - 16). 

To combine the component distributions for the composite 
distribution pertinent to Figure (73), we take the mean value of 
the distributions in the overlapped regions. It is believed that with 
the satellite speed ratio which is in the order 10, the suggested means 



Figure 73 
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of treating the mutual interference effect should be accurate within ±10%. 
(b) charged particle distribution - To evaluate the ion and electron 
distribution and the electric potential around the satellite, Figure (.73), 
by the use of Figures (19 - 58), we must first determine the following 
parameters : 

(1) R/d where d, the Debye length is given in Table I, d^sO. 5 cm, 
hence, the ratios of R/d. for components A, B, and C are estimated 
to be in the order R/d?^20 if R = 10 cm. Since the results will 
not be sensitive to the values of R/d as long as R/d 1, hence the 
order R/d = 20 is close enough. 

(2) The temperature ratio T /T. must be estimated for practical 

estimation, take T /T. = 1, unless T /T. is known to be much 

e l e i 

larger than 1. There is, however,, only a limited case of calculations 
given for T / T. / 1. See e.g. Figures(59 - 7 2). 

(3) To determine the speed ratio V / v. , we use Table I to obtain 

the ion thermal speed v ^ v^(neutral) ~ 1.5 x 10 5 , then to estimate 

V /v. from the given value for V , the satellite speed, 
o 1 o 

(4) The surface potential (<j)g) , as well as the other electric potential 
<f>£ given in the figures are all normalized by the quantity kT./ e. 

Knowing T from Table I (usually it is assumed T = T , the kinetic 
temperature of neutral particles), the quantity kT./e can be calculated, 
e.g. with T = 2000° K, kT./e = 1,72. (cj>* = eF/kT.) 

1 O o 1 

With the four parameters (R/d , T^/T. , V/w , cj>g) prescribed 
as shown above, we can use Figures(l,9 - 58) to determine the ion 
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density and potential behind the components A, B } ~ and C, 

respectively. -The corresponding electron density n can be calculated 

* <})❖ 

from the use of Maxwell- Boltzmann law, namely n = n e , where 

e o 

n^ is the ambient electron density which is shown in Table I. 

To combine the component solutions for the composite configuration 
of Figure(73) we again take the mean value of corresponding quantities 
in the overlapped regions to reduce the mutual interference effects. 

It is believed that the results so obtained should be accurate within 
±15% of the true solution with precisely prescribed parameters of 
R/d , Tg/T^ , V/v^j and <$>g for the actual configuration. 

The plasma sheath for the components A, B, and C can be 
determined by first segmenting the surfaces into small inclined planes 
and then using the Figures(17) and (18) for the ion density (n.) and 

electric potential (6*). The electron density can be determined again 

, 

by n = n e T . 
e o 
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Legend of Figures 19-71 

The cases of n^- and ^-distribution in the wake of a spher< 

or a cylinder presented in Figures 19-71 are specified by four 

parameters (R/d, T /T. } V /v., ed> /kT.) and a siibscript "s" or 

e i o l S i 

"c" to denote "sphere" or "cylinder", respectively. The distributions 
are given in two different coordinate systems: 

i) downstreams sections[see Figure l] : z/R = 1, 2, 3, 4 with 

(etfj/kT^) vs (p/R) or (n^/n^) vs (p/R) for sphere; 

(e4>/kTJ vs (y/R) or (n./n^) vs (y/R) for cylinder. 

ii) radial directions: 9 = 1°, 27°, 46°, 60°[ see Figure i] with 

A - i_ 

(e0/kT.) vs [(p 2 + z 2 ) 2 /R] or {n./n ) vs E(p 2 + z 2 ) 2 /r] 
i i o 

for sphere; 0 = 0°, 27°, 45°, 60° with (eejs/kT.). vs [ (y 2 + z 2 ) 2 /r] 

i 

or (n./nj vs [ (y 2 + z 2 ) 2 /R] for cylinder. 



Figure I - Sphere or Cylinder 
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Appendices B and C (Self-consistent solutions of charged particle- 

and potential-distribution in the near wake of a satellite). 

The work reported in Appendix C ("Near Wake of the Rarefied 
Plasma Flows at Mesothermal Speeds, " A. I. A. A. Preprint No. 
68-169) is an improved version of that in Appendix B ("On the 
Electro dynamic Wake in a Rarefied Plasma, " Rarefied Gas Dynamics,' 
5th. Symposium, vol. 2). The results presented in Figures (19 - 71) 
are based on a further improved version which will be described in 
detail in a forthcoming Ph.D. thesis (H. Jew, University of Michigan, 
1968). The essential improvement over that in Appendix C includes 
a more rigorous expression for the third constant of motion (I 3 ) in 
the case of wake behind a sphere and the removal of a computational 
error in the case of wake behind a cylinder. 



